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How robust are the constraints on cosmology and galaxy 
evolution from the lens-redshift test? 
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ABSTRACT 

The redshift distribution of galaxy lenses in known gravitational lens systems provides 
a powerful test that can potentially discriminate amongst cosmological models. How- 
ever, applications of this elegant test have been curtailed by two factors: our ignorance 
of how galaxies evolve with redshift, and the absence of methods to deal with the effect 
of incomplete information in lensing systems. In this paper, we investigate both issues 
in detail. We explore how to extract the properties of evolving galaxies, assuming that 
the cosmology is well determined by other techniques. We propose a new nested Monte 
Carlo method to quantify the effects of incomplete data. We apply the lens-redshift 
test to an improved sample of seventy lens systems derived from recent observations, 
primarily from the SDSS, SLACS and the CLASS surveys. We find that the limiting 
factor in applying the lens-redshift test derives from poor statistics, including incom- 
plete information samples, and biased sampling. Many lenses that uniformly sample 
the underlying true image separation distribution will be needed to use this test as 
a complementary method to measure the value of the cosmological constant or the 
properties of evolving galaxies. Planned future surveys by missions like the SNAP 
satellite or LSST are likely to usher in a new era for strong lensing studies that utilize 
this test. With expected catalogues of thousands of new strong lenses, the lens-redshift 
test could offer a powerful tool to probe cosmology as well as galaxy evolution. 

Key words: gravitational lensing - galaxies: evolution ~ galaxies: luminosity func- 
tion, mass function - cosmological parameters - cosmology: observations - cosmology: 
theory 



1 INTRODUCTION 

Gravitational lensing statistics have now been used to map 
the mass distribution in galaxies (Blandford & Narayan 
1992; Narayan & Bartelmann 1996; Kochanek 2004) as well 
as to constrain cosmological parameters (Cheng & Krauss 
1998; Chae 2003; Maoz 2005). Since the discovery of the 
first multiply imaged quasar (Walsh, Carswell & Weymann 
1979), well over a hundred such systems have now been dis- 
covered in various wave-bands, ranging from the optical to 
the radio. This progress is attributed to several dedicated 
on-going all sky surveys like the CLASS (Myers et al. 2003) , 
the SDSS (York et al. 2000), the Two-Degree Field Galaxy 
Redshift Survey (2dFGRS; CoUess et al. 2001) and more re- 
cently the SLACS (Bolton et al. 2006). Consequently, there 
has also been significant progress in analytical and statisti- 
cal studies of lenses. Many sophisticated methods are now 
available to model strong lensing systems using both para- 



* e-mail: capelo@astro.yale.edu 
f e-mail: priya@astro.yale.edu 



metric and non-parametric mass distributions (e.g. Rusin, 
Kochanek & Keeton 2003; Saha & Williams 2003). 

The concept of optical depth to lensing (ODTL) was 
proposed to study strong lensing statistics by Turner, Os- 
triker & Gott (hereafter TOG, 1984): they presented an ana- 
lytic calculation of the lensing probability of distant quasars 
by intervening galaxy lenses and the role of selection ef- 
fects therein. Since the lensing probability depends on the 
comoving volume element, the ODTL test can be used to 
constrain cosmological parameters by comparing the num- 
ber of expected lenses to the number of observed ones. Using 
the ODTL test, Kochanek (hereafter K96, 1996) obtained 
limits on the cosmological constant from the statistics of 
gravitational lenses using a number of completed quasar sur- 
veys (e.g. Snapshot Survey, the ESO/Liege survey, the NOT 
survey, the HST GTO survey, the FKS survey), lens data, 
and a range of lens models. The formal limit obtained was 
JIa < 0.66 at 95% confidence in fiat cosmologies, which in- 
cluded the statistical uncertainties in the number of lenses, 
galaxies, quasars, and the parameters relating galaxy lumi- 
nosities to dynamical variables. This value is in contrast to 
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what is now well established by WMAP observations of the 
CMB (e.g. Spergel et al. 2006) and high redshift SN la obser- 
vations (Riess et al. 1998; Perlmutter et al. 1999). These ob- 
servations have in fact led to what is currently referred to as 
a 'cosmic concordance' model (Ostriker & Steinhardt 1995) 
- the A-CDM model (with h ~ 0.7, Oa ~ 0.7, D.m 0.3 
and Qk — 0.0) - as the most widely accepted description of 
the Universe. It has been argued by K96 that their retrieved 
low value for SIa could be due to dust obscuration in a large 
fraction of lensing galaxies: however, a hundred times more 
dust is needed to change the expected number of lenses by 
a factor of two. Given this extreme value, dust is clearly 
not the dominant source of systematic errors. By tabulating 
various sources of error and the limitations imposed on the 
accuracy of the determination of Qa, K96 speculated that 
the assumptions on the velocity dispersion function of lenses 
might be a significant source of error. Reviewing previous 
estimates of the cosmological constant derived from strong 
lensing statistics Maoz (2005) concludes that the discrepan- 
cies might be due to possibly a lower lensing cross section for 
ellipticals galaxies than assumed in the past. Maoz (2005) 
argues that the current agreement between recent model cal- 
culations and the results of radio lens surveys may be for- 
tuitous, and due to a cancellation between the errors in the 
input parameters for the lens population and the cosmology, 
as well as input parameters for the source populations. 

In the quest to determine the correct underlying cos- 
mological model by placing better and tighter constraints 
on f^A, strong gravitational lensing has not been the most 
reliable technique. Systematic errors have plagued the lens- 
ing analysis, leading to contradictory results for the derived 
values of the cosmological constant in a flat Universe (see, 
for example, Maoz & Rix 1993, K96, Chae et al. 2002). 
These contradictory results were primarily caused by: small 
number statistics due to the shortage of observed lens sys- 
tems; assumptions about the relationship between luminosi- 
ties and masses of galaxies; scatter in the empirical relation 
between mass and light; and observational biases, mainly 
the magnification bia^. An explicit relation between mass 
and light is required for the lensing analysis in the absence 
of independent mass estimates for the lensing galaxies. The 
luminosity of galaxies is converted into a mass distribution 
(which is the relevant quantity to model lensing effects) us- 
ing a density profile, which is parametrized via the velocity 
dispersion. Statistics of strong lenses and any cosmological 
constraints thereby obtained depend on the assumed veloc- 
ity dispersion function (VDF) of galaxies. 

Kochanek (hereafter K92, 1992), devised a test, the 
'lens-redshift test', which circumvented the magnification 
bias since it does not involve computing the total ODTL. 
This test relies on the computation of the differential op- 
tical depth to lensing with respect to the angular critical 
radius r. The probability distributions of lens redshifts z 
with a given angular critical radius, [{dr / dz) / T]{r) , are eval- 
uated. However, this quantity still required knowledge of the 
VDF of lensing galaxies, which was inferred (hence IVDF) 
by combining the Schechter luminosity function with an em- 
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pirical relation between luminosity and velocity dispersion, 
the Faber-Jackson and the TuUy-Fisher relations for early- 
type and late-type galaxies, respectively. The lens-redshift 
test depends on cosmological parameters as well as on galaxy 
evolution parameters. Therefore, it can be used to constrain 
the former by fixing the latter, or vice versa. With the as- 
sumption of no evolution, K92 derived Qa < 0.9. 

Ofek, Rix & Maoz (hereafter ORM, 2003) revived the 
lens-redshift test (K92). They applied it to a larger sample 
of lens systems than were available to the K92 analysis, us- 
ing the CLASS and SDSS surveys. Their study also included 
a re-derivation and generalisation of the lens-redshift test, 
which incorporated mass and number density evolution of 
lens galaxies. They explicitly included the redshift evolution 
of the characteristic velocity dispersion and evolution of the 
number density of galaxies. The limit obtained by ORM for 
a flat Universe, assuming no mass evolution of early-type 
galaxies between z = 0-1, was Qa < 0.95 at the 99% con- 
fidence limit. Turning things around, and fixing the cosmo- 
logical model to S7a ~ 0.7 and Qm = 0.3, they determined 
galaxy evolution parameters and found d log j^^a (z) /dz — 
-0.10lo;e andd logj^^nt^iz) /dz = -1-0. 7li;2, where a* and n^, 
are the characteristic velocity dispersion and number density 
of lensing galaxies, respectively. 

Mitchell et al. (hereafter MKFS, 2005) focused instead 
on the ODTL test (TOG). In addition to using a larger 
sample than the one used by K96, they included the evo- 
lution of the VDF in amplitude and shape, based on the- 
oretical galaxy formation models, and used the measured 
velocity distribution function (MVDF) for early-types from 
the SDSS (Sheth et al. 2003). MKFS found ^a = 0.74-0.78 
for a flat Universe prior and a limit S7a < 0.86 at the 95% 
confldence limit. Including the effects of galaxy evolution, 
they found Qa = 0.72-0.78 and a limit Qa < 0.89 at the 
95% confldence limit. 

The consequence of using the MVDF versus the IVDF 
in the determination of Qa is one of the key questions we 
address in this work. The IVDF and MVDF diflfer at high 
luminosities/ velocity dispersions (see Fig. Al in the Ap- 
pendix). The scatter of the Faber-Jackson relation was a 
predominant source of uncertainty in the previous studies 
(cf. ORM), leading to a systematic underestimation of the 
number of objects with large velocity dispersions. 

In this paper, we investigate the lens-redshift test in de- 
tail and re-examine the uncertainties that limit its use as a 
powerful discriminant between cosmological models, as well 
as its potential to constrain galaxy evolution models. We ap- 
ply this to a new enlarged sample of lenses. This is done for 
the flrst time using the measured velocity dispersion func- 
tion from SDSS although we compare and reproduce the 
results of ORM using the inferred velocity dispersion func- 
tion. In addition, we consider the effect of incomplete lensing 
information on the retrieval of cosmological parameters with 
a new nested Monte Carlo method. 

The outline of the paper is as follows. In section 2 we de- 
flne the lens-redshift test and compare the use of the IVDF 
and the MVDF on the determination of both Qa and galaxy 
evolution parameters. In section 3 we describe the new ex- 
panded sample and, in section 4, we present the results of 
the application of the lens-redshift test to our sample. We 
present a new Monte Carlo method to quantify the effect of 
incomplete lensing information in section 5, by constructing 
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realizations of several biased subsamples. We conclude with 
a discussion of our results and their implication for future 
observational surveys. 



2 THE LENS-REDSHIFT TEST FORMALISM 

2.1 Methodology using the inferred velocity 
dispersion function 

We follow the notation introduced by TOG and K92 in defin- 
ing the optical depth to lensing and the lens-redshift test, 
respectively. The differential optical depth to lensing per 
unit redshift is the differential probability dr that a line of 
sight intersects a lens at redshift z in traversing the path dz 
from a population of lensing galaxies with comoving num- 
ber density hl- Mathematically, for a source this is simply 
the ratio of the differential light travel distance cdt to its 
mean free path between successive encounters with galaxies 
I/ulS, 



dr 
dz 



= nL{z)S 



Eiz)il + z)' 



(1) 



where the comoving number density of lensing galaxies is 
given by wl = n*(l + 2:)''; n* is the average number density 
of lensing galaxies; S is the cross section for multiple imaging 
of a background point source; th ~ c/Hq is the Hubble 
radius and E{z) = {flM{l + z)'^ + ^^(1 + z)^ + Qa)^ . The 
cross section for multiple imaging S is given hy S — tw^Dl ^. 
We initially assume n*, the characteristic luminosity L* and 
the characteristic velocity dispersion cr* of lensing galaxies 
to be constant with redshift, although we will later relax this 
assumption and allow for time evolution. 

Our analysis is restricted to early-type and SO galaxies 
as lenses and it is assumed that they can be modelled as 
singular isothermal spheres (SIS) 0. With the assumptions 
stated above, we can write the angular critical radius r as: 
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where Dls and Ds are the angular diameter distances be- 
tween the lens and the source and between the observer and 
the source, respectively, and Je is a parameter that takes 
into account the difference between the velocity dispersion 
of the mass distribution Om and the observed stellar veloc- 
ity dispersion a: Om ~ fE(^- Modelling galaxies as singular 
isothermal spheres, the characteristic central velocity dis- 
persions (which are typically unmeasured for most lenses), 
are drawn from the VDF. 

In this subsection, we relate the luminosity distribution 
to the Faber- Jackson law and construct the IVDF. Using 
the Schechter function fit to model the luminosity function 
of lensing galaxies, 

dn n* / L \ " r i 1 
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and the Faber- Jackson relation, L/L* = {a/a-^y , to relate 
the luminosity to a velocity dispersion, combining these two 
equations we derive the IVDF: 
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Combining the IVDF with equation (1), the differential 
optical depth can be written as 
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dz {da/a-i. 
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Defining r* = 47r(a*/c)^ and using 

r 1 Ds 
/b^ Dls 



gives us the IVDF lens-redshift test equation. 
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where and = IQ-k'"^ Ui^r [a / c)^ are constants. 

Incidentally, we note that this is slightly different in 
form from K92, as we compute dr /dzdr, whereas K92 com- 
pute dr I dz(dr I r . Both calculations then proceed to nor- 
malise with respect to r; this gives identical results only 
when a single population of galaxies is considered. The value 
of r* depends on ct*, which in turn varies as a function of 
the morphological type considered. We include both ellipti- 
cals and SO galaxies in our analysis, whereas K92 considered 
only ellipticals. 

To obtain constraints on galaxy evolution, we consider 
the following scaling relations where 

n*(^)=naO^^ (8) 

L4^) = Lao'?^ (9) 

a,(2) =craO^", (10) 

n*, L* and cr^ are the characteristic values at zero redshift, 
and P, Q and U axe constants (as in equations (9), (10) and 
(11) in ORM). Incorporating these, the IVDF lens-redshift 
test equation (7) becomes 

'Dls 1 
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Note that Q, the evolution parameter of the luminosity 
in equation (9), does not appear in the equation above. 



^ A SIS has a mass distribution given by p(r) = /2-KGr^ , where 
CT is constant with radius r. This density profile is a very good fit 
for elliptical and SO lensing galaxies. In fact non-singular isother- 
mal spheres and truncated isothermal spheres give very similar 
fits to lensing data (e.g. Rusin et al. 2003). 



2.2 Methodology using the measured velocity 
dispersion function 

In this section we rewrite the differential optical depth as 
a function of the measured velocity dispersion function for 
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early-type galaxies from observations circumventing the use 
of the Schechter luminosity function and the Faber-Jackson 
relation. We now use the functional form (fitted to SDSS 
observations) of the MVDF taken from Sheth et al. (2003) 
and rewrite it in a form that makes for easy comparison with 
the IVDF (equation (4) above): 
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exp 
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Substituting and simplifying the expression for the 
differential optical depth as done previously, we obtain the 
MVDF lens-redshift test equation, 
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where = W'K^nJrH^io'l/c)'^ and = 47r(o-i/c)^ are 
constants. 

Once again, we consider the case where the parameters 
ni, Z/t and ai for the MVDF evolve with redshift, following 
equations similar to (8), (9) and (10), and on substituting 
we find 
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dzdr 
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exp 
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There is no reason to assume that cr* and cr* (or n* and 
n*/) evolve differently, therefore, we set U = U' (or P = P'). 
We note here that cr* is just a parameter in the IVDF and 
the MVDF fits and due to the different functional forms 
for the parametrizations of the velocity dispersion function, 
their values for the IVDF and MVDF can be and are in fact, 
found to be quite different. 



3 DEFINING THE NEW LENS GALAXY 
SAMPLE 

Following ORM, our first lens sample is primarily drawn 
from the CASTLES (Muiioz et al. 1998) data baseEl which, at 
present, contains 82 class 'A' (certain), 10 class 'B' (likely), 
and 8 class 'C (dubious) gravitational lenses, making for a 
total sample size of 100 systems. We ignore the 13 class 'B' 
binary quasars from the CASTLES lists, eight of which have 
image separations greater than 4" and would be discarded as 
likely being cluster-assisted rather than due to field galaxies. 
In the remaining 5 lenses nearby group galaxies are impli- 
cated in determining the separations, therefore we discard 
them as well. 

To get a handle on potential biases in the sample we 
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have grouped the systems by their discovery technique into 
three categories: targeted optical discoveries, which were tar- 
gets selected based on the lensed source optical emission 
and includes such surveys as the first HST snapshot survey 
(Maoz et al. 1993) and quasars selected from the Calan- 
Tololo (Maza et al. 1996), Hamburg-ESO (Wisotzki et al. 
2000) and SDSS quasar surveys (York et al. 2000); tar- 
geted radio discoveries, which were targets selected based 
on the lensed source's radio properties and includes the 
JVAS/CLASS (Myers et al. 2003), PMN (Winn et al. 2002), 
and MG (Bennett et al. 1986) lens searches; and miscel- 
laneous discoveries for systems discovered either serendip- 
itously (such as the HST parallel field discoveries; Rat- 
natunga, Griffiths & Ostrander 1999) or based on system 
properties other than that of the lensed source and typically 
that of the lensing galaxy. 

Lenses discovered based on the properties of the source 
ought not to harbour a bias in terms of the redshift of 
the lensing galaxy, although they suffer from magnification 
bias. However, systems discovered because of the lens or 
surrounding environment will naturally favour low-redshift 
lenses. All systems in the 'Miscellaneous Discoveries' fall un- 
der this category, which includes systems discovered because 
of the properties of the lensing galaxy such as Q2237-(-030 
(Huchra et al. 1985) and CFRS03.1077 (Crampton et al. 
2002), systems discovered based on properties of the lens- 
ing galaxy's environment such as RXJ0921-I-4529 (Mufioz et 
al. 1998) , and systems discovered serendipitously from HST 
pointings such as the HST Medium Deep Survey lensing 
candidates (Ratnatunga, Griffiths & Ostrander 1999) and 
HDFS2232509-603243 (Barkana, Blandford & Hogg 1999), 
the latter of which are characterised by deflector emission 
that is either comparable to or dominates over the back- 
ground source. 

We also exclude systems inappropriate for our lens 
model of isolated and elliptical/SO lensing galaxies. These 
include systems with multiple lensing galaxies of comparable 
luminosities (and therefore likely comparable halo masses) 
such as HE0230-2130 (Wisotzki et al. 1998), B1359-H54 
(Myers et al. 1999) and B2114+022 (Augusto et al. 2001). 
We also exclude cluster-assisted systems such as Q0957-f 561 
(Young et al. 1980). 

Although we are ignoring entire surveys, this ought not 
to introduce biases in the sample. These various cuts de- 
tailed above leave a total of 42 systems in our sample 
Al detailed in Table Al with complete redshift information 
(source and lens redshifts). We have estimated the size of the 
deflector's critical radius for the remaining systems using a 
simple Singular Isothermal Sphere model (SIS) plus external 
shear using the gravlens software of Keeton (2001). Rela- 
tive image positions with respect to the lensing galaxy were 
obtained from either the CASTLES compilation or from the 
reference in column (11) of Table Al in the Appendix. For 
double systems, we use the reddest flux ratio measured be- 
tween lensed components (typically either HST/F160W or 
the radio flux ratio if the system is radio-loud) as the fifth 
constraint required by the model, which ought to minimize 
flux ratio contamination from microlensing-induced variabil- 
ity. For systems with ring morphology, the critical radius 
of the lens was obtained from the corresponding model re- 
ported in the cited reference. 

Finally, we add the SLACS lenses (Bolton et al. 2006) 



How robust is the lens-redshift test? 5 



6 - 




2 4 6 

critical radius (arcseconds) 



Figure 1. Histogram of tiie critical radii of lenses for the various 
samples considered in this work: sample Al (brown and orange); 
Bolton et al. (2006) SLAGS sample of 28 lenses (19 confirmed 
lenses and 9 candidates) (blue). Note that the SLAGS lenses have 
systematically larger critical radii as a consequence of their spec- 
troscopic selection from the SDSS. 



to construct our largest sample comprising seventy systems. 
The SLAGS survey has proven to be a very efficient HST 
Snapshot imaging survey for new galaxy-scale strong lenses. 
The targeted lens candidates were selected by Bolton et al. 
(2006) from the SDSS database of galaxy spectra for hav- 
ing multiple nebular emission lines at a redshift significantly 
higher than that of the target SDSS galaxy. The survey is op- 
timized to detect bright early-type lens galaxies with faint 
lensed sources. The key advantage of this selection tech- 
nique is that it provides a homogeneously selected sample 
of bright early-type lens galaxies. However, given the size of 
the fibres in SDSS this sample is biased toward large sep- 
arations (the consequences of this bias are discussed later) 
compared to other surveys. This is clearly seen in the his- 
togram plotted in Fig. 1. All the 28 SLAGS lenses found to 
date are tabulated in Table A2 of the Appendix. Of these 
lenses 19 are confirmed with multiple images and the re- 
maining 9 lenses are candidates. We note here that both 
the confirmed and unconfirmed candidates from the SLAGS 
are included in our calculations. Although the SLAGS is a 
biased sample, we include all 28 lenses in our analysis to 
illustrate the effect of better statistics given the success of 
the adopted selection strategy. The SLAGS selection favours 
gravitational lenses that have typically larger Einstein radii 
(as seen clearly in the histogram of image separations in 
Fig. 1) by virtue of the selection of spectroscopic candi- 
dates for imaging follow-up. For the purposes of the current 
analysis it clearly strengthens the results of this work, i.e. 
biased sampling of the image separation distribution pro- 
vides biased values for galaxy evolution parameters. Once 
the survey has finished, the selection function will be very 
well determined, which will enable more careful use of this 
sample for galaxy evolution studies. This strategy has been 
extremely successful so we feel compelled to showcase this 
sample. SLAGS E/SO lenses appear to be a random sub- 
sample of the luminous red galaxies sample of the SDSS, 
only skewed toward the brighter and higher surface bright- 
ness systems. While the environments of some of the lenses 
are complicated by the existence of nearby galaxies, by and 



large it is a 'clean' sample where the image separations are 
determined primarily by a single elliptical/SO lens. Includ- 
ing the 28 SLAGS lenses to the sample Al gives us a final 
tally of 70 lenses, all with complete information, that defines 
our sample A2. In Fig. 1, we plot the distribution of crit- 
ical radii for observed lenses in sample Al and the SLAGS 
lenses, which together constitute our sample A2. Since sev- 
eral samples will be used in the paper, we list them here 
for clarity: Sample Al - our updated version of the 
ORM sample I, with a total of 42 systems; Sample 
A2 - our new, enlarged sample that contains Sample 
Al and 28 new SLACS lenses; Sample B - our mock 
sample of a 100 lenses with complete information; 
and Sample C: our truncated sample Al with 10% 
of the largest separation lenses removed. 



4 ANALYSIS AND RESULTS 

In this section, we compare the discriminating power of the 
MVDF versus the IVDF, in constraining both cosmological 
and galaxy evolution parameters. We study the recovery bias 
in the extraction of (i) cosmological constraints with U = 
P = for the various compiled lens samples, as well as of 
(ii) galaxy evolution parameters, by fixing the cosmological 
parameters. Finally, we assess the impact of incompleteness 
of lens data in the recovery of the cosmological constant. 

As noted in section 2, we normalise the lens-redshift 
probability distribution with respect to the optical depth r. 
Therefore, parameters which appear simply as multiplying 
constants in the distributions do not impact our comparison. 
Such parameters include the Hubble radius vh = c/Ho and 
the average number density of lensing galaxies n* (IVDF) 
and nj (MVDF). Our comparison of the IVDF and MVDF 
is not affected by the value of /e- This parameter relates 
the velocity dispersion of the dark matter to that of the 
stars. TOG set it to (3/2)^^^, other studies (e.g. Narayan 
& Bartelmann 1999) suggested using values smaller than 1. 
Recent results from the SLAGS survey suggest that fs ^ 1, 
i.e. the lens model velocity dispersions are fairly close to 
the measured stellar velocity dispersion within an effective 
radius (Treu et al. 2006). Therefore, we take /_b = 1 in 
this work. The default cosmological model is taken to be a 
Friedmann-Robertson- Walker {QA + i^M+^K = 1) A-GDM 
flat Universe, with flu ~ 0.3, Hk ~ 0.0 and Qa ~ 0.7. 

The parameters that affect our analysis are the ones 
defining the VDFs: ct*, 7, a, a'^, a', 13'. For the IVDF: 
a is the faint-end slope in the Schechter luminosity func- 
tion, and 7 is the Faber-Jackson power-law index. For the 
MVDF: a' is the low-velocity power-law index, and is 
the high-velocity exponential cutoff index of the distribu- 
tion. Following ORM, the valued (ji^e, ?^*so, c"*, Q, 7) = 
(0.0039/1^ Mpc-^ 0.0061/1^ Mpc-^ 225 km s-\ -0.54, 4) are 
used for the IVDF. Following MKFS fit to the MVDF, we 



^ CTi, = 225 km s~^ is the characteristic velocity dispersion for 
elliptical galaxies only; SO galaxies have a characteristic velocity 
dispersion of 206 km . When varying Oi, in our calculations, 
we also vary the characteristic velocity dispersion of SO galaxies 
as (T*x (206/225). 
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taks0 (n*',crl,Q',/3') = (0.0041/^3Mpc-^ 88.8 km s-\ 6.5, 
1.93). 

Note that late-type galaxies can in principle also be in- 
corporated into the IVDF analysis easily, by simply replac- 
ing the Faber- Jackson relation with the TuUy-Fisher rela- 
tion. We can substitute the Faber- Jackson exponent 7 and 
characteristic cr* with the corresponding TuUy-Fisher rela- 
tion parameters. Although late-type galaxies are more nu- 
merous than early-type galaxies, they tend to have lower 
masses and therefore do not contribute significantly to the 
total optical depth to lensing. Due to the strong dependence 
of the lensing cross section on the velocity dispersion, this 
causes late-type galaxies in general to be inefficient lenses. 
Besides, late-type galaxies are not included in the determi- 
nation of the MVDF from SDSS data. So in this work, for 
consistent comparisons we restrict ourselves to elliptical and 
SO lenses. 

We compute the differential optical depth {dT/dz)/T for 
each individual lens with measured separation and known 
source redshift for our two samples: sample Al (42 lenses) 
and sample A2 (70 lenses). We then determine the proba- 
bility distribution of the redshift of the lens using equation 
(7) (IVDF) and equation (13) (MVDF), given the observed 
image separation; the measured source redshift; a given cos- 
mology and galaxy evolution model (in this case assuming 
no evolution: U = P = Q = 0, h = 0.7, D.m = 0.3, Qk = 0.0 
and Q.A = 0.7). If the choice of underlying cosmological pa- 
rameters, primarily Qa in this case, corresponds to the true 
value, the peak of the probability distribution Zp ought to 
be close to the measured lens redshift z;, in a good num- 
ber of cases. These lens redshift probability distributions are 
shown in Fig. A2 (in the Appendix) for all the lenses in our 
sample A2. The plot illustrates that the MVDF and IVDF 
yield near identical probability distributions for the lens red- 
shifts. Choosing different values of Ha shifts these inferred 
probability distributions: this is illustrated in Fig. 2 for one 
lens (B0218+357), for values of = 0.2,0.4,0.7 and 1.0 
(keeping ^Ik = 0.0). However, we notice a systematic ef- 
fect: Zp the peak redshift of the probability distribution is 
skewed slightly lower for the MVDF compared to the IVDF 
for almost the entire sample A2. This can be qualitatively ex- 
plained as arising due to the different asymptotic behaviours 
of the MVDF and the IVDF (see Fig. Al in the Appendix). 

4.1 The maximum likelihood method 

Now, we use the samples to statistically play the game in 
both directions: (i) constrain the geometry of the Universe 
with galaxy evolution parameters fixed, and (ii) to constrain 
the galaxy evolution model with cosmological parameters as 
knowns. 

We use a maximum likelihood estimator in our sta- 
tistical analysis of lens redshift distributions. The lens- 
redshift test equations give the probability distribution of 
lens redshifts as P{zi\{X}, Zs,r), normalised to unity, where 
{X} = {Qa , U, P} is the set of cosmological and galaxy evo- 
lution model parameters, and {zs,r) are source redshift and 

^ We note that Sheth et al. (2003) choose a higher value for n*'; 
moreover, n^.^ + ^*S0 7^ 'i*': but the choice of these parameters 
does not affect our analysis. 
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Figure 2. Dependence of the lens-redshift distribution on cos- 
mology. We show the dependence on the value of Ha explicitly 
for the lens B0218-I-357 using the MVDF (thick lines) and IVDF 
(thin lines) respectively, for a range of values: Qa = 0.2 (black, 
dotted), Qa = 0-4 (red, short-dashed), Qa = 0.7 (green, solid) 
and Qa = 1.0 (blue, long-dashed). The vertical dashed line marks 
the position of the observed lens redshift zi = 0.69. The peak Zp 
of the distribution increases as Qa increases from 0.2 to 1.0, due 
to the increase in the cosmological comoving volume. 

lens angular critical radius priors for a given system. The 
likelihood estimator L for the entire sample of A'^ systems is 
then: 

N 

Li{X})=l[p(zi\{X},Zs,r). (15) 

i = l 

The quantity L is computed to quantify the consistency 
of all measured lens redshifts for the entire ensemble of lens 
systems for any given geometry and galaxy evolution model. 

We then compute the maximum value of L, fixing the 
galaxy evolution parameters to obtain constraints on the 
cosmology {{X} — {Qa}), and then fixing the cosmology to 
obtain constraints on galaxy evolution {{X} — {U,P}). 

As pointed out by ORM, the lens-redshift test is more 
sensitive to the galaxy mass evolution parameter U com- 
pared to the galaxy number evolution parameter P. This 
can be understood by considering the limit when P ~ 0: a 
negative U decreases the most probable value for the lens 
redshift and narrows the probability distribution. In con- 
trast, the number evolution parameter only affects the peak 
value but does not affect the overall shape of the probability 
distribution. 

4.1.1 Constraints on cosmology 

We proceed to obtain constraints on JIa, keeping galaxy 
evolution parameters fixed, using the Friedmann-Robertson- 
Walker cosmology, and imposing Qa + Qm + Qk ~ 1 with 
Qk ~ 0.0. Unless otherwise stated, we assume U — P = Q = 
0, corresponding to the case of no evolution in the galaxy 
population either in mass or number. The lens-redshift test 
equations (7) and (13) are used in this instance and the 
likelihood as described above is constructed and maximized. 

A projection of the likelihood surface along the Qa axis 
for sample Al is shown in Fig. 3. In the upper panel the like- 
lihood function is calculated using the IVDF for sample Al. 
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Several values of ct* are shown for completeness. Assuming 
a value of ct* — 225 km s~^ for elliptical galaxies (ORM), 
we obtain the following limits on the cosmological constant: 
Qa = 0.55to'2o If'" confidence. This value is determined by 
taking the median as the central value and "ruling out" the 
leftmost ~ 16% and rightmost ~ 16% of the total integral 
('median' method). Alternatively, if we take the mode as the 
central value and determine the threshold value of the likeli- 
hood for which the integral under it comprises ~ 68% of the 
total integral ('mode' method), we obtain a slightly higher 
value of = Omto ll; this method is similar to what 
ORM did, except that they assumed a normal distribution 
('normal' method). Their assumption turns out to be quite 
reasonable as we recover the same resuhs {Qa = 0.60to ll) 
applying their method to our sample Al. For the IVDF, we 
note that the error bars we obtain are much smaller than 
in ORM: this is certainly due to the larger number of lens 
systems in our sample Al (the ORM sample I had 15 lenses 
compared to 42 in our sample Al). Although our error bars 
are smaller, we note that the numbers quoted here are for 
a single value of ct* — 225 km s~^ and these determinations 
of r^A using the 'mode', 'median' and 'normal' methods are 
entirely consistent with each other within the errors. 

The corresponding results for sample Al using the 
MVDF are also shown in Fig. 3 (lower panel): in this case 
as well, several values of ai are plotted to present the trend 
clearly. Assuming a'^ = 88.8 km s~^ (Sheth et al. 2003, 
MKFS), we obtain a value for the cosmological constant of 
= 0.62lo;i7 using the 'median' method, Qa = 0.67lo:ig 
using the 'mode' method, and Q,a = 0.671q;i5 using the 
'normal' method. Again, we note that these quoted values 
are for a single value of = 88.8 km s~^ and once again 
the constraints on Qa using these three different criteria are 
consistent with each other given the errors. Even with the 
improvement of using the MVDF compared to earlier work, 
the sensitivity to Qa in the lens-redshift test is low, as seen 
clearly by the fact that using the ±lcr range on recov- 
ers values of JIa varying from 0.0 to nearly 1.0 - the full 
available range. 

Unsurprisingly, the recovery of Qa using the MVDF and 
the IVDF yields very similar values, as these functions differ 
only at the extremely high (a* > 380 km s~^) velocity dis- 
persion tail as shown in Fig. Al in the Appendix. A marked 
difference between the IVDF and MVDF shows up in the 
velocity range of 380-400 km s~^, which is characteristic 
of cD galaxies. Strong lensing events from such galaxies are 
difficult to model, as their position at the centre of clusters 
causes the events to be assisted by additional smoothly dis- 
tributed dark matter in their vicinity. Since we excluded all 
such systems in our sample Al, it is not surprising that the 
inferred value of the cosmological constant using the IVDF 
and MVDF are in good agreement. 

We plot the corresponding results for the larger sample 
A2 in Fig. 4 employing the MVDF. Once again we plot the 
projection of the likelihood for various values of a'^. For a'^ = 
88.8 km s~^, we now find Qa = 0.86lQ Qg for the 'median' 
method; Q,a = 0.87to'o6 for the 'mode' method and Qa = 
0-87to:o6 for the 'normal' method. 

We find that the recovered value of Qa is higher from 
the sample A2 (the mode value is shifted by about 0.25). 
Sample A2 does include a higher proportion of larger separa- 
tion lenses (clearly seen in Fig. 1). This indicates a potential 




"a 

Figure 3. The projection of the likelihood surface along the Qa 
axis is shown here for several values of the characteristic velocity 
dispersion cr* (IVDF; upper panel) and tr^ (MVDF; lower panel) 
for sample Al. The numbers denote the characteristic velocity 
dispersion in km s"'^, increasing in value from the top down in 
both panels as indicated by the arrows. For tr* = 225 km s~^ 
(black thick line, upper panel), we obtain Qa = 0-60^q :[g. For 
(jj^ = 88.8 km s~^ (black thick line, lower panel), we obtain 
Qj^ = The likelihood curves are very shallow and con- 

sequently the error bars are rather large. The MVDF and IVDF 
constraints on Qa are in good agreement. 

systematic bias that skews recovery of Qa, that is sensitive 
to how well the 'true' separation distribution is sampled. To 
obtain robust constraints on Qa with the lens-redshift test 
not only do we need large samples but we also need lenses 
that accurately reflect the true underlying distribution of 
image separations. We note here that the SLAGS lenses are 
included to clearly demonstrate this bias as their image sep- 
arations are skewed to larger values as a consequence of the 
selection technique. 

Four key results emerge from these plots: first, the lens- 
redshift test is not very robust in constraining the value of 
the cosmological constant with current samples. This was al- 
ready suggested by K92, but we demonstrate it more clearly 
here even with two notable improvements: a larger sample 
of lenses and the use of the MVDF. The likelihood curve 
is very shallow and consequently the error bars are rather 
large. Second, the MVDF and IVDF results are compara- 
ble, therefore the inefEcacy of the lens-redshift test does not 
appear to stem from systematics arising from the use of the 
IVDF. Third, is the notable sensitivity of constraints on Qa 
to the parameter cr, . The value of cr, emerges in the fit of a 
functional form to the observed velocity dispersion function 
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and depends on the completeness of the measurement, i.e. 
adequate samphng of the high and low velocity dispersion 
tail for observed galaxies. Finally, inclusion of the SLAGS 
lenses (19 confirmed lenses + 9 candidates) with relatively 
larger separations pushes the recovered Ha to higher values. 
The finite number of lens systems is clearly implicated here 
as evidenced in the error bars on ^Ia and is a key limitation. 
In conclusion, as we show in the next section, while a large 
number of lens systems will go a ways toward increasing the 
robustness of this test in the future, it is crucial to simulta- 
neously sample the separation distribution uniformly. 

4-1-2 Constraints on galaxy evolution 

We now investigate galaxy evolution using the lens-redshift 
test, with n*, L*, cr*, n^, and (T'^, varying with redshift 
according to equations (8), (9) and (10) and their primed 
versions. The equations used are (11) and (14). Fixing the 
cosmological model io h = 0.7, Q-m = 0.3, Q-k ~ 0.0 and 
Qa = 0.7, we determine U and P- As outlined in Section 
4.1, the likelihood function is constructed fixing SIa and then 
maximizing to obtain constraints on U and P. 

We obtain constraints on U and P for various samples: 
the ORM sample I lenses, sample Al and sample A2 all of 
which are plotted in Fig. 5. We calculate the U-P contours 
for the ORM sample I, applying our analysis methods and 
using the MVDF. The MVDF was not available at the time 
of the ORM analysis. Our calculation of the U and P pa- 
rameters using the IVDF is in very good agreement with 
their results. The orientation and calibration of the confi- 
dence level contours agree. Although the contours of all the 
samples overlap quite well, the difference in the peak values 
of U and P determined for our samples Al, A2 and ORM 
sample I is significant. 

The likelihood results in the U-P plane for sample Al, 
using the IVDF and the MVDF, are shown in the upper and 
the lower panels of Fig. 5, respectively. We obtain a maxi- 
mum in L at {[/ = 0.11, P = -1.40} for the IVDF and at 
{U = 0.10, P = -1.24} for the MVDF, again showing no 
significant dependence on the choice of velocity dispersion 
function employed. However, we note that the contours close 
along the (7-axis for the MVDF case compared to the IVDF. 
Therefore, using the IVDF to calculate the likelihood low- 
ers the sensitivity to mass evolution. For our sample A2 the 
maximum value for L is at {U = 0.32, P = —1.60} for the 
IVDF and at {U = 0.28, P = -1.57} for the MVDF. For 
the ORM sample I: the maximum value for L is found to lie 
at {U = -0.08, P = 0.44} for the IVDF and at {U = -0.07, 
P = 0.85} for the MVDF. The likelihood peak moves toward 
increasingly positive values of U for sample A2 compared to 
Al and the ORM sample I. This indicates a strong sensitiv- 
ity to the fraction of large separation lenses. Sample A2 has 
a larger proportion of those and therefore predicts stronger 
mass evolution for the lens ensemble. 

The fact that U and P have opposite signs is consistent 
with mass conservation: we have either a larger number of 
lower mass galaxies (P > 0, C/ < 0) or fewer more massive 
galaxies {P < 0, U > 0) in the past compared to today. 
However, the case (P < 0, C/ > 0) is in conflict with the 
currently accepted hierarchical model of galaxy formation 
with bottom-up assembly of structure. 

Our primary conclusions on deriving galaxy evolution 




Figure 4. The projection of the likelihood surface along the Qa 
axis is shown here for several values of the characteristic veloc- 
ity dispersion cr'^ for the full sample A2 using the MVDF. The 
numbers denote the characteristic velocity dispersion in km s~^, 
increasing in value from the top down as indicated by the ar- 
rows. For For a'^ = 88.8 km s''^ (black thick line), we obtain 
Qa = O.ST^QQg. The likelihood curves are very shallow and con- 
sequently the error bars are rather large. The sample A2 clearly 
yields values of Qa that are systematically higher than sample 
Al. 



parameters are: (i) we reproduce the trends reported by 
ORM for their sample I when we use the IVDF as they 
did; (ii) we find the likelihood peak position to be insensi- 
tive to the choice of IVDF vs MVDF for the ORM sample 
I; (iii) we recover U and P values consistent within the er- 
rors for all our samples; (iv) with a larger number of lenses 
(as in sample A2) we obtain slightly increased sensitivity 
to P compared to ORM sample I; (v) the likelihood peak 
shifts systematically to higher U values for sample A2 which 
contains a higher proportion of large separation lenses com- 
pared to the sample Al. To summarise, there is a notable 
observational bias in recovering mass evolution that depends 
strongly on how well the underlying true separation distri- 
bution is sampled in detected lenses. 



4-1-3 Investigation of systematic observational biases 

Unbiased lens surveys are needed to sample uniformly the 
full distribution of separations (Kochanek 1993) in order to 
apply the lens-redshift test to constrain galaxy evolution 
parameters as found above. If the sample is slightly skewed 
toward larger separations, biased values of the galaxy evolu- 
tion parameters are retrieved. To investigate this issue fur- 
ther, we create a mock sample of a hundred lenses (sample 
B), all with complete information (i.e. source redshift, lens 
redshift, and image separation) known assuming no evolu- 
tion, to try and understand the observational biases that 
likely affect our analysis. We randomly assign the source 
redshift from a normal distribution of redshifts centred at 
z = 2 and with a dispersion of 1. We then randomly assign 
the angular critical radius, as half of the lens separation, 
from the probability distribution of image separations given 
by Kochanek (1993), where we set 7 = 4 and a = — 1 in 
their equation (4.10) for a fiat Universe. We calculate the 
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Figure 5. Comparison of constraints on galaxy evolution for different samples: the derived U-P contours for sample Al (blue, solid 
lines), sample A2 (black, dot-dashed lines) and the ORM sample I (red, dashed lines). In the upper panel, the IVDF was used for the 
calculation and in the lower panel the MVDF. The three contours shown for each sample are the Itr, 2(T and 3cr confidence levels and the 
cross marks the no-evolution locus {U = 0, P = 0). While there is some degree of overlap for the various samples it is clear that the peak 
value of U - the mass evolution parameter - shifts toward more positive values consistent with sample A2 having a larger proportion of 
more massive lenses, in good agreement with the fact that sample A2 does have a higher fraction of large separation lenses. The error 
bars on the ORM sample I are larger since it has only 15 systems whereas sample Al has 42 systems and sample A2 70 systems. We 
note that within the errors the values of U and P obtained for different samples are in good agreement. 



differential optical depth distribution for each lens, using 
equation (7), and randomly pick a lens redshift from it. 

For the full mock catalogue (sample B), the maximum 
value for the likelihood L is obtained at {U — —0.03, P = 
—0.57}. The input {U — 0, P — 0) parameters values are 
not exactly recovered due to finite-sample variance. Several 
subsamples of sample B were then evaluated, after cutting 
the sample based on source redshifts and image separations. 
Creating a subsample of 90 lenses discarding the 10 highest 
source redshift systems, we find that the maximum value of 
the fikelihood L shifts to {U = -0.09, P = -0.21}. Then 
for a subsample of 90 lenses generated by discarding the 10 
lowest source redshift systems we find that the peak is at 
{U = -0.03, P = -0.54}. 

To examine the additional sensitivity to the number 
of lenses as well, we now construct a subsample from sam- 
ple B of 58 lenses with the lowest source redshifts and find 
{U = —0.07, P — —0.80} and for the subsample of 58 lenses 
with the highest source redshifts the evolution parameters 
are found to be {U = -0.02, P = -0.48}. These resuhs are 
shown in the lower panel of Fig. 6: no significant system- 
atic bias is introduced on selection by source redshift. The 



only difference seen in these subsamples is the effect of the 
variation in the number of systems: the contours are more 
extended for the 2 cases with 58 systems compared to the 
cases with 90 systems. 

We then cull sample B on the basis of image separations 
and the results for these subsamples are shown in the upper 
panel of Fig. 6. The key systematic that we study in further 
detail is the role of biased sampling of the image separa- 
tion distribution. The mock catalogue generated above was 
now cut based on lens angular critical radius r. Once again 
biased subsamples were generated to preferentially sample 
larger and smaller separation systems. First, we constructed 
a subsample of 90 systems discarding 10 smallest separation 
systems: for this instance the maximum value of the likeli- 
hood lies at {U = -0.03, P = -0.49}. Picking now a further 
90 systems discarding the 10 largest separation systems, we 
find a different maximum at {U = —0.02, P = —0.62}. 
Similarly, making a more extreme selection, we pick 58 sys- 
tems from the mock discarding 42 of the largest separa- 
tion systems. This is our extreme biased sample skewed to 
small separations. For this subsample we find {U — —0.14, 
P — —0.19}. Finally, for a mock with 58 of the largest sep- 
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Figure 6. Effect of biased lens samples and the recovery of U and P illustrated for the IVDF case. Upper panel [Biased image separations]: 
U-P contours for our full mock sample B and several subsamples selected by cuts in image separations (i) lowest 90% separations (blue, 
dot-dashed thin lines), (ii) lowest 58% separations (red, dotted thin lines), (iii) largest 90% separations (blue, solid thin lines), (iv) largest 
58% separations (orange, dashed thin lines). Lower panel [Biased source redshifts]: U-P contours for subsamples now cut on the basis of 
source redshifts (i) lowest 90% Zs (blue, dot-dashed thin lines), (ii) lowest 58% Zs (red, dotted thin lines), (iii) highest 90% Zs (blue, solid 
thin lines), (iv) highest 58% Zs (orange, dashed thin lines). The full mock sample B results (black, solid thick lines) are shown in both 
panels for comparison. The three contours shown for each sample are the la, 2a and 3(7 contour likelihood lines and the cross marks the 
location of the no-evolution locus (C/ = 0, P = 0). A clear systematic bias is introduced on selection by lens separation. In particular, 
the velocity dispersion evolution parameter dramatically shifts from near-zero values (£/ ~ 0) to negative values {U ~ —0.1), when the 
largest separation lenses are removed from the sample. 



aration lenses (this constitutes our extreme biased sample 
toward large separations) we find {U = —0.10, P = 1.92}. 
Our analysis clearly indicates the presence of a systematic 
bias introduced on selection by lens separation. This effect 
is especially seen clearly in the smaller subsamples (with 58 
systems). In particular, the velocity dispersion evolution pa- 
rameter dramatically shifts from near-zero values {U ~ 0) 
to negative values {U ~ —0.1), when the highest separation 
lenses are removed from the sample. We see clearly from 
Fig. 6 that lens data comprising biased sampling of the un- 
derlying image separation distribution introduce a system- 
atic shift in the recovered values of the galaxy evolution 
parameters, whereas data with biased source redshifts yield 
unbiased estimates of U and P. 

Galaxy evolution parameters are thus extremely sensi- 
tive to observational biases in the separation distribution of 
lens systems. For ground based optical surveys and high res- 
olution HST surveys there are optimal separations that are 
detected. Lens systems found in ground based surveys are 



likely skewed to larger separations than those found in HST 
surveys. 

Having narrowed the plausible source of the systematic 
bias, we re-do the analysis making a similar cut on our ob- 
served lens sample Al to verify our finding. We remove the 
five (10%) largest image separation lenses thus creating the 
sample C. We now compare the recovery of U and P for 
this biased sample with samples Al and A2, sample B and 
the ORM sample I. The maximum value of L is obtained 
for sample C at {[/ = -0.03, P = -0.51} (IVDF) and at 
{U = -0.01, P = -0.40} (MVDF). In Table 1 we list the 
positions of the peak values of the likelihood function in the 
U-P plane, for our samples Al, A2, B, ORM sample I and 
sample C, for the IVDF and MVDF. The full results are 
shown in Fig. 7. The resultant trend clearly demonstrates 
the strong bias now replicated with cuts in the data intro- 
duced by artificially removing large image separation sys- 
tems. Our result that incompleteness in image separations 
is a serious current limitation in using strong lensing statis- 
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Figure 7. Illustrating the systematic bias arising from poorly sampling the true image separation distribution: the U-P contours for 
ORM sample I (red, dotted lines) and our samples Al (blue, solid lines), A2 (black, dot-dashed lines), C (orange, dashed lines), and B 
(black, solid lines) are shown for the IVDF (upper panel) and MVDF (lower panel). The three contours shown for each sample are the 
Icr, 2(T and 3iT confidence levels, with the cross marking {U = 0, P = 0) denoting the no-evolution case. The systematic bias induced in 
our sample A2 arises from a larger fraction of large image separation lenses. 



Sample U=d log]^Q(T*(2:)/c(2 P=d logiQnt(z)/dz 



A2 (IVDF) +0.32 -1.60 

Al (IVDF) -1-0.11 -1.40 

B (IVDF) -0.03 -0.57 

C (IVDF) -0.03 -0.51 

ORM I (IVDF) -0.08 4-0.44 



A2 (MVDF) +0.28 -1.57 

Al (MVDF) +0.10 -1.24 

B (MVDF) +0.00 -0.71 

C (MVDF) -0.01 -0.40 

ORM I (MVDF) -0.07 +0.85 



Table 1. Positions of the peak values of the likelihood function 
in the U-P plane, for samples A2, Al, B, C, and ORM I, IVDF 
and MVDF, in order of U-peak position. When comparing sam- 
ples A2, Al and C, we confirm the same trend we observed with 
the mock subsamples: the velocity dispersion evolution parameter 
shifts to less positive values as we remove the highest separation 
lenses. 



tics has also been pointed out by Oguri (2005) and Oguri, 
Keeton & Dalai (2005). 



5 THE EFFECT OF INCOMPLETE LENS 
DATA ON THE RETRIEVAL OF 
COSMOLOGICAL PARAMETERS 

Previous lens-redshift test analyses have differed on how to 
handle systems with incomplete redshift information. K96 
included an estimate of the probability of failing to measure 
a system's lens redshift for systems lacking such a measure- 
ment. ORM take a more pragmatic approach by discarding 
all systems with Za > 2.1, arguing that systems below that 
redshift are mostly complete. The former approach is made 
difRcult by the many variables that can prevent a success- 
ful redshift measurement (surface brightness of the lensing 
galaxy, galaxy contrast with respect to the magnified source 
images, observing conditions during an actual measurement 
attempt), while the latter approach ignores higher-redshift 
systems that do have complete redshift information. Such 
systems are likely to show the strongest sensitivity to cos- 
mological or evolution effects that are sought after in the 
first place. 

The approach we adopt here is to marginalise over 
systems with incomplete redshift information using nested 
Monte Carlo simulations. Let Nc be the number of systems 
with complete redshift information and Nu be the number 
of systems with unmeasured lensing redshifts. For a given 
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parameter set {X}, we can assign lens redshifts for the Nu 
sample by drawing from P{zi\{X}, Zs, r) which gives a sam- 
ple of lens redshifts With {2;,^} fixed, we obtain the 
absolute likelihood L{{X}) for the combined Nc + Nu sam- 
ple. The procedure is then repeated Nmc times with each 
iteration using a different set of {^z.u}. This gives an aver- 
age absolute likelihood < L{{X}) > and a corresponding 
scatter 5L{{X}) for the given set of model parameters {X}. 
The scatter in the absolute likelihood estimate shrinks to 
zero as Nu — > 0, and can be interpreted as a measure of the 
uncertainty in the absolute likelihood because of the incom- 
plete sample. The entire procedure can then be repeated for 
a different set of model parameters. 

We argue that this is an attractive method for several 
reasons. First, it does not ignore existing redshift informa- 
tion for any system, either within the complete or the incom- 
plete sample. This results in as large a sample size as possible 
and helps to minimize small-number effects that are tradi- 
tionally present in lensing statistics. Second, the question of 
handling biases present in the incomplete sample is made 
objectively by marginalising over the entire sample rather 
than imposing an artificial cut on, say, the source redshifts. 
And third, it allows one to quantify the effects that the in- 
complete sample has on the accuracy of likelihood analysis 
through SL{{X}). This last point can be used to explore how 
the precision of the model parameters can be measured by 
future changes in either the complete or incomplete sample 
size. 

We performed a nested Monte Carlo simulation of our 
sample Al {Nc = 42), adding a set of ten mock lens sys- 
tems with known image separations, source redshifts and 
unknown lens redshift (A^„ = 10), to make up a total of 
52 lens systems. We fixed all galaxy evolution parameters 
{U = P = Q = 0, a* = 225 km s~^) and varied cosmo- 
logical parameters, after setting S7a + + ^Ik ~ 1 with 
Qk = 0.0. Therefore, the parameter set was taken to be 
{X} = {JIa}. a projection of the (un- normalised) likelihood 
surface along the Qa axis is shown in Fig. 8 for the IVDF 
case: we show the results for Qa values steps of of 0.1 only, 
and for Nmc = 100. The blue line is the result for the sample 
Al. For each fixed value of 57a, each of the points represents 
the value of the likelihood for one of the Nmc = 100 different 
possible sets of the combined Nc 4- Nu sample. Approximat- 
ing the likelihood distributions with Gaussian functions of 
mean \i =< L({X}) > and dispersion a — SL{{X}), the 
average absolute likelihood < L{{X}) > for the combined 
Nc + Nu sample follows the pattern given by the complete 
system (blue line), being almost flat between Oa ~ 0.0 and 
Qa ~ 0.7. Therefore adding even a small number of systems 
with incomplete information reduces further the sensitivity 
to QaE 

The effect of incomplete information is to further di- 

^ The likelihood distributions have very large dispersions indicat- 
ing the lack of robustness in determined values of the cosmological 
constant. Below we enumerate some typical values of the mean 
H =< L{{X}) > and dispersion a = SL{{X}) in the likelihood 
stepping through a grid of Ha values that illustrates this point: 
= 0.0: A* ^ -12.6, cr ~ 7.9; Qa = 0.1: ~ -10.0, a ~ 4.7; 
Qa = 0.2: fj. ~ -11.2, cr ~ 7.4; Qa = 0.3: ^ -11.8, a ~ 8.0; 
Qa = 0.4: ^ ~ -11.5, cr ~ 8.6; Qa = 0.5: A* ^ -10.2, cr ~ 7.8; 
Qa = 0.6: A* ^ -9.7, ct ~ 8.1; Qa = 0.7: fi ~ -13.3, cr ~ 9.9; 
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Figure 8. The projection of the (un-normalised) likelihood sur- 
face along the Qa axis is shown here for hundred independent 
realizations of 52 lens systems (IVDF case), comprising 42 sys- 
tems from sample Al and 10 mock systems. All galaxy evolution 
parameters being fixed {U = P = Q = 0, cr^ = 225 km s~^), each 
dot represents the value of the likelihood L for a particular value 
of Qa and for a particular set of 52 lenses. The blue line is the 
result for sample Al (42 lenses) only. For each value of Qa, the 
distribution of 100 likelihood values can be approximated by a 
Gaussian function with a mean fj, =< L{{X}) > and a dispersion 
a = SLiiX}). 

lute the efficacy of constraints on cosmological parameters. 
While we have argued here that the lens-redshift test with a 
small lens sample with complete information is insufficient, 
we further find even with a small number of systems in a 
large sample with incomplete information, we lose sensitiv- 
ity to Qa- 



6 CONCLUSIONS AND DISCUSSION 

We investigate the lens-redshift test to assess its robustness 
in constraining cosmological and galaxy evolution param- 
eters. We apply the test to a much improved lens sample 
compared to earlier work by K92 and ORM. Moreover, we 
also use the observationally determined velocity dispersion 
function (MVDF), instead of relying on the IVDF. MKFS 
also used the MVDF, but they applied it to the ODTL test, 
which is more affected by observational biases - mainly the 
magnification bias - than the lens-redshift test considered 
here. Finally, we develop a new nested Monte Carlo analysis 
to quantify the effects of incompleteness on the accuracy of 
retrieving S^a- 

Our results suggest that the lens-redshift test is not 
particularly robust in the determination of either cosmo- 
logical parameters or galaxy evolution parameters with the 
currently available samples. We conclude this after careful 
analysis of 70 lens systems and generating several mock cat- 
alogues. First, we fix galaxy evolution parameters to con- 
strain Oa: in this instance very weak constraints are ob- 
tained. Moreover, despite using the MVDF for the first time 



Qa = 0.8: fi ~ -12.2, cr ~ 5.5; Qa = 0.9: fj. ~ -18.9, cr ~ 5.3; 
Qa = 1.0: fj. ~ -35.3, cr ~ 6.9. 
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in this test our results do not differ significantly from ear- 
lier work. When we do the converse, i.e. fix the cosmology 
and look for constraints on galaxy evolution, we find that 
the results on the evolution parameters are too sensitive to 
the choice of sample, implying a very strong dependence 
on the observational bias introduced by lens separations. 
Finally, the limit of the precision with which the value of 
the cosmological constant can be determined due to lack of 
complete information has been assessed. We find that even 
a small number of systems with incomplete information in 
a large sample can further reduce the significance of the 
already weak constraints on the cosmological constant. In 
fact, systems with incomplete information add more noise 
than signal. For the purposes of constraining cosmological 
parameters incomplete-redshift information systems are best 
excluded. With the small number of systems available at the 
present time, such a strong cut is not feasible, however with 
the expected large number of new lenses from future sur- 
veys the statistics will permit stricter selection of optimal 
systems. 

The lens-redshift test is clearly affected by lens selec- 
tion effects. An obvious observational strategy for the fu- 
ture would be to observe hundreds of new lenses, that fairly 
sample the full distribution of separations. Such samples are 
expected from the large area surveys to be performed by fu- 
ture instruments like SNAP and the LSST. These large sam- 
ples with hundreds/thousands of lenses at several redshifts 
would allow us to better quantify the lens sample selection 
bias. Moreover, ideally lenses for use in the lens-redshift test 
need to be relatively "clean" , that is, they should not belong 
to groups, where presence of additional deffectors/nearby 
galaxies could affect the image separation and therefore pro- 
vide skewed lens image separation distributions that will 
in turn bias results. It is becoming increasingly clear from 
the study of individual lens environments that there are 
clearly no truly isolated lenses. However, what is important 
from the point of view of the lens-redshift test is that the 
neighbouring perturbers are not massive enough to signif- 
icantly alter the image separations to within observational 
positional accuracies. Obviously these accuracies depend on 
whether space based data or ground based data is available 
for new lens systems. From large proposed future surveys 
which all involve deep imaging, lens systems with many per- 
turbers need to be culled. It appears from simulations that 
systems that likely require of the order of 10 - 20% exter- 
nal shear are still viable candidates for the lens-redshift test 
(Oguri, Keeton & Dalai 2005) . 

Future lens samples should ideally include both early 
and late type galaxies and span a large redshift range, in 
order to constrain galaxy evolution parameters. Finally, all 
lens systems to be useful should have complete information, 
since even a small fraction of incomplete systems would sig- 
nificantly decrease the efficacy of constraints on parameters. 
In the near future, hundreds/thousands of new lens systems 
will be discovered by these upcoming new instruments. Si- 
multaneously with many on-going and planned ambitious 
surveys to study galaxy evolution, progress is likely to come 
from better knowledge of galaxy evolution. The lens-redshift 
test, which is currently unable to give decent constraints on 
cosmology and galaxy evolution because of poor statistics, 
could eventually prove to be a very profitable means to con- 



strain cosmological parameters and galaxy evolution models 
robustly. 
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Table Al. The columns are: (1) lens number; (2) lens name; (3) R.A. [h:m:s], ,12000.0; (4) Dec. [d:m:s], J2000.0; (5) source redshift, 2^; 
(6) lens redshift, zj; (7) critical radius, r, in arcseconds; (8) grade for the likelihood that the object is a lens: A=rd bet my life, B=rd 
bet your life, and C=rd bet your life and you should worry (CASTLES); (9) Number of images corresponding to each source component, 
E means extended and R means there is an Einstein ring (CASTLES); (10) sample: Al- Sample Al in this paper, C- Sample C in this 
paper, I- Sample I in ORM, O- Targeted optical discoveries, R- Targeted radio discoveries, M- Miscellaneous discoveries; (11) references. 
The numbers shown are the ones actually used in the computations. When computing ORM Sample I, we used the numbers from ORM 
Table Al (calculating the critical radius, r, as half of the separation). 
List of references: 

1 - Wisotzki et al. (2004); 2 - Ofek et al. (2006); 3 - Rusin et al. (2003); 4 - Morgan et al. (1999); 5 - Wisotzki et al. (2002); fi - Hagen 
& Reimers (2000); 7 - Johnston et al.(2003); 8 - Kneib, Cohen & Hjorth (2000); 9 - Kochanek et al. (2000); 10 - Lubin et al. (2000); 11 
- Eigenbrod et al. (2006a); 12 - Schechter et al. (1998); 13 - Lehar et al. (2000); 14 - Lacy et al. (2002); 15 - Hewett et al. (1994); 16 - 
Surdej et al. (1997); 17 - Eisenhardt et al. (1996); 18 - Lidman et al. (2000); 19 - Weymann et al. (1980); 20 - Tonry (1998); 21 - Sluse 
et al. (2003); 22 - Oguri et al. (2005); 23 - Eigenbrod et al. (2006b); 24 - Burud et al. (2002a); 25 - Wisotzki et al. (1996); 26 - Burud 
et al. (2002b); 27 - Cohen, Lawrence & Blandford (2003); 28 - Wiklind & Combes (1995); 29 - Lawrence et al. (1995); 30 - Tonry & 
Kochanek (1999); 31 - Fassnacht & Cohen (1998); 32 - Myers et al. (1999); 33 - Patnaik et al. (1992); 34 - Lehar et al. (1993); 35 - Treu 
& Koopmans (2003); 36 - Winn et al. (2002); 37 - Langston et al. (1989); 38 - Wikhnd & Combes (1996); 39 - Sykes et al. (1998); 40 - 
Warren et al. (1996); 41 - Ratnatunga et al. (1998); 42 - Gregg et al. (2000); 43 - Fassnacht et al. (1996); 44 - Morgan et al. (2001); 45 - 
Eigenbrod et al. (2007); 46 - Surdej et al. (1993); 47 - Morgan et al. (2005). 
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Number 


Name 




Za 




<Ta 


Status 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 




SDSS 


lOO'^T^'i 21 — nQ422n 1 


0.632 


0.195 


265 + 10 

^\JtJ _1_ 


Q 


2 


SDSS 




0.524 


0.332 


QQO + 23 


Q 


3 


SDSS 




0.581 


0.322 


310 ± 15 


Q 


4 


SDSS 




0.446 


0.194 


231 ± 16 


uc 


5 


SDSS 


TfiQi 2(15 sn+nn2Qni i 


0.324 


0.164 


313 ± 12 


Q 


g 


SDSS 


TnQ532n 424-52n54S 7 


0.467 


0.131 


207 ± 14 


UC 


7 


SDSS 




0.470 


0.241 


299 ± 16 


Q 


g 


SDSS 


T0Q'SQ44 07-1-041017 


0.535 


0.126 


212 ± 12 


Q 


9 


SDSS 


1102551 '51—0(13517 4 


0.276 


0.159 


247 ± 11 


Q 


10 


SDSS 


Till 73Q fiO+05341 3 Q 


0.823 


0.229 


279 ± 21 


UC 


11 


SDSS 


1120540 434-4Q102Q 3 


0.481 


0.215 


iiOO _1_ xu 


uc 


12 


SDSS 


1125028 25-I-05234Q 


0.795 


0.232 


254 ± 14 


Q 


13 


SDSS 


11 9^1 '^'^ 70 — 090^^0'^ 1 


0.784 


0.224 


216 + 23 


Q 


14 


SDSS 




0.449 


0.233 


263 ± 17 


uc 


15 


SDSS 


1133045 53 — 014841 6 


0.712 


0.081 


178 ± 09 


Q 


16 


SDSS 


Tl 40998 91 -l-fi'^91 '^'^ 


0.481 


0.205 


275 + 15 


Q 


17 


SDSS 


114201 5 85+601 Ql 4 8 


0.535 


0.063 


194 ± 05 


Q 


18 


SDSS 


J154731.22+572000.0 


0.396 


0.188 


243 ± 11 


UC 


19 


SDSS 


J161843.10+435327.4 


0.666 


0.199 


257 ± 25 


C 


20 


SDSS 


J162746.44-005357.5 


0.524 


0.208 


275 ± 12 


C 


21 


SDSS 


J163028.15+452036.2 


0.793 


0.248 


260 ± 16 


C 


22 


SDSS 


J163602.61+470729.5 


0.675 


0.228 


221 ± 15 


uc 


23 


SDSS 


J170216. 76+332044.7 


0.436 


0.178 


239 ± 14 


uc 


24 


SDSS 


J171837.39+642452.2 


0.737 


0.090 


270 ± 16 


c 


25 


SDSS 


J230053. 14+002237.9 


0.464 


0.229 


283 ± 18 


c 


26 


SDSS 


J230321. 72+142217.9 


0.517 


0.155 


260 ± 15 


c 


27 


SDSS 


J232120.93-093910.2 


0.532 


0.082 


236 ± 07 


c 


28 


SDSS 


J234728.08-000521.2 


0.715 


0.417 


330 ± 50 


uc 



Table A2. List of all lenses belonging to Bolton et al. (2006) sample. The columns arc: (1) lens number; (2) lens name; (3) source 
redsliift, Zs; (4) lens redsliift, Zf, (5) velocity dispersion measured within an aperture Ca, in km s^"'" and (6) the status of the lens (C 
- confirmed; UC - unconfirmed). The numbers shown arc the ones actually used in the computations (only the central value of (To was 
used). The critical radius r was then computed using equation (2) and h = 0.7, Qa = 0.7, VIm = 0.3, and Qk = 0.0. Note that the 19 
confirmed lenses and 9 candidates are denoted with a C and U C respectively in the final column. 




Figure Al. The measured velocity dispersion function (MVDF) and the inferred velocity dispersion function (IVDF) for early-type 
galaxies. The functional forms plotted here are derived from fits provided in equation (23) of Mitchell et al. (2005). The solid curve is 
the MVDF and the dashed curve is the IVDF. 
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Figure A2. The lens-redshift distribution {{dT/dz)/T vs z) for all 70 lenses from Sample Al (numbered in the same order as presented 
in Table Al in the Appendix) and the Bolton at al. (2006) SLACS sample (numbered in the same order as presented in Table A2 in the 
Appendix): calculated using the IVDF (red, dot-dashed line) and the MVDF (blue, solid line); with the galaxy evolution and cosmological 
parameters set to U = P = Q = 0, h = 0.7, Q/^ = 0.7, Qj^l = 0.3, and Qj^ = O-O- The vertical dashed line marks the position of the 
observed lens redshift zi. The peak Zp of the probability distribution is skewed to slightly higher zi for the IVDF compared to the MVDF, 
in most lenses. 



